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Abstract. We first study the properties of the Fuchsian ordinary differential 
equations for the three and four-particle contributions x' 3 ' an d of the square 
lattice Ising model susceptibility. An analysis of some mathematical properties 
of these Fuchsian differential equations is sketched. For instance, we study 
the factorization properties of the corresponding linear differential operators, 
and consider the singularities of the three and four-particle contributions x' 3 ' 
and x' 4 ^ i versus the singularities of the associated Fuchsian ordinary differential 
equations, which actually exhibit new "Landau-like" singularities. We sketch the 
analysis of the corresponding differential Galois groups. In particular we provide 
a simple, but efficient, method to calculate the so-called "connection matrices" 
(between two neighboring singularities) and deduce the singular behaviors of 
X^ 3 ' and x' 4 ' ■ We provide a set of comments and speculations on the Fuchsian 
ordinary differential equations associated with the n-particle contributions X 
and address the problem of the apparent discrepancy between such a holonomic 
approach and some scaling results deduced from a Painleve oriented approach. 
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1. Introduction 

It is a great pleasure to write this paper, in honour of the sixtieth birthday of Prof. 
A.J. Guttmann, on the mathematical structures deduced from series expansion of the 



| Dedicated to A.J. Guttmann on the occasion of his 60-th birthday. 



Beyond series expansions. 



2 



susceptibility of the Ising model. One of us (JMM) has profited from many years of 
fruitful collaboration and correspondence with him and his group, and many times 
from his friendly hospitality. 

1.1. Notation and a few results on the Ising susceptibility. 

Since the (monumental) work of T.T. Wu, B.M. McCoy C.A. Tracy and E. Barouch 
[1], it has been known that the expansion in n-particle contributions to the zero field 
susceptibility of the square lattice Ising model at temperature T can be written as an 
infinite sum : 

oo 

x{T)=Y J X in) {T) (1) 

71=1 

of (n — l)-dimensional integrals [2, 3, 4, 5, 6, 7]. This sum is restricted to odd 
(respectively even) n for the high (respectively low) temperature case. While the 
first contribution in the sum, x^\ is obtained directly without integration, and the 
second one, x^ 2 \ is given in terms of elliptic integrals, no closed forms for the higher 
order contributions were known, despite the well-defined forms of the integrands in 
these (n — l)-dimcnsional integrals. It is only recently that the differential equations 
for an d have been found [8, 9, 10, 11, 12]. In the following we will use, as 
in [8, 9, 10, 11], the following normalization : 

x in) = (l-O . £(„) (2) 

in order to focus on the X^'s which arc only functions of the variable w = 
s/(l + * 2 )/2. 

As far as singular points are concerned (physical or non-physical singularities in 
the complex plane), and besides the known s = ±1 and s = ±i singularities, B. Nickel 
showed [6] that ^( 2n+1 ) is singular! for the following finite values of s = sh(2J/kT) 
lying on the \s\ =1 unit circle (to = k = excluded): 

( s ) u k u m 

u 2n+1 =1, -n < to, k < n (3) 

When n increases, the singularities of the higher-particle components of x( s ) 
accumulate on the unit circle \s\ = 1. The existence of such a natural boundary 
for the total susceptibility x( s )> shows that x( s ) ls n °t D-finite (not holonomic) as a 
function§ of s. 

A significant amount of work had already been performed to generate isotropic 
series coefficients for x*™' (by B. Nickel [6, 7] up to order 116, then to order 257 by 
A.J. Guttmannff et al). More recently, W. Orrick et al. [15, 16], have generated 
coefficients^ of the susceptibility of the Ising model [15, 16] x( s ) U P to order 323 (resp. 
646) for high (resp. low) temperature series in s, using a quadratic double recursion 
(non-linear Painleve difference equations) for the correlation functions [15, 16, 17, 

J The singularities being logarithmic branch points of order e 2n (™+ 1 )- 1 . l n {e) with e = 1 — s/s, 
where s; is one of the solutions of (3). 

§ See in particular, the paper of I. Enting and A.J. Guttmann [13]. One has probably the same 
situation for the three-dimensional Ising model [14]. 
ff Private communication. 

J The short-distance terms were shown to have the form (T — T c ) p ■ (log\T — T c \) q with p > q 2 . 
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18, 19]. As a consequence of this remarkable non-linear Painleve difference equation 
and its associated double recursion, the computer algorithm had a 0(N 6 ) polynomial 
growth for the calculation of the series expansion instead of the exponential growth 
that one might expect at first sight. 

However, in such a non-linear, non-holonomic, Painleve- oriented approach, one 
obtains results directly for the total susceptibility x(s) which do not satisfy any finite 
order linear differential equation. It thus prevents the easy disentanglement of the 
contributions of the various holonomic x^'s. 

On the other hand, the individual X s do satisfy finite order linear differential 
equations, leaving some hope to understand the susceptibility from a deeper knowledge 
of the mathematical structures of the successive differential equations of the n-particle 
sequence of the x^'s. 

2. The Fuchsian ODE's for x (3) and x (4) - 
2.1. The Fuchsian ODE for x (3) . 

With an original methodfj [9] that allows us to write yf® as IUU Y integrated sums, 
we generated a long series of 490 terms. Using a dedicated program we searched for 
the finite order linear differential equation with polynomial coefficients in the variable 
w = s/2/(l + s 2 ), by steadily increasing the order. We finally succeeded [8] in finding 
the following linear differential equation of order seven for x' 3 ^ with only 359 terms 
of our long series: 

E °» • - ( 4 ) 

with : 

a 7 = 

w 7 ■ (1 - w) (l + 2u;) (l-4io) 5 {l + Awf (1 + 3w + 4w 2 ) ■ P 7 {w), 
a 6 = w 6 -(l-Awfil + iwf -P 6 {w), (5) 
a 5 = w 5 -(l-Awf (1 + 4w) ■ P 5 (w), 

a 4 = • (1 -Aw) 2 ■ P 4 (w), a 3 = w 3 ■ (1 - Aw) ■ P 3 (w) 

a 2 = w 2 ■ P 2 {w), ai = w ■ P\{w), a = Pq{w) 

where [8] : 

P 7 (w) = 1568 + 15638 w - 565286 w 2 - 276893 w 3 + ••• 

+ 158329674399744 w 27 + 39582418599936 w 28 (6) 

The polynomials P 7 (w), P&(w) ■ ■ ■, Po(w) are polynomials of degree respectively 
28, 34, 36, 38, 39, 40, 40 and 36 in w, and are given in [8]. The singular points of this 
differential equation correspond to the roots of the polynomial corresponding to the 
highest order derivative F^(w). 

The linear differential equation (4) is an equation of the Fuchsian type since there 
are no singular points, finite or infinite, other than regular singular points. With 

(j It is worth noting that in the details of this method we encountered many times the occurrence 
of generalizations of hypcrgcomctric functions to several variables, extremely similar to the ones 
occurring in Feynman diagrams [20] (Appel, Lauricella, Kampe de Feriet, ...). We will discuss these 
quite interesting questions elsewhere. 
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w = 


s = 




P = 


9,3,2,2,1,1,1 


w = -1/4 


s = -1 




P = 


3,2,1,0,0,0,-1/2 


to = 1/4 


s = 1 




P = 


1,0,0,0,-1,-1,-3/2 


w = -1/2 


1 + s + s 2 = 




P = 


5,4,3,3,2,1,0 


w = 1 


2 - s + 2s 2 = 




P = 


5,4,3,3,2,1,0 


l + 3to + 4to 2 = 


(2s 2 + s+ l)(s 2 + s - 


h2) = 


P = 


5,4,3,2,1,1,0 


1/to = 


1 + s 2 = 




P = 


3,2,1,1,1,0,0 


to = top, 28 roots 


s = sp, 56 roots 




P = 


7,5,4,3,2,1,0 



Table 1: Critical exponents for each regular singular point, wp is any of the 28 roots 
of Pyiw). We have also shown the corresponding roots in the s variable§. 

this property, using the Frobenius method [21, 22], it is straightforward to obtain the 
critical exponents, from the indicial equations, in to, for each regular singular point. 
These are given in Table 1. 

We have actually found two remarkable rational and algebraic solutions of (4), 
namely : 

o 

s = w s = w (7) 

1 l-4to' 2 (l-4to)Vl-16w 2 

We also found a solution behaving like to 3 , that we denote S3 : 

5*3 = to 3 + 3w 4 + 22to 5 + 74to 6 + 417tt; 7 + 1465to 8 

+ 7479 w 9 + 26839 w 10 + ■ ■ ■ (8) 
and associated with the physical solution '■ 

S 9 = to 9 + 36 w 11 + 4 w 12 + 884 to 13 + 196 to 14 

+ 18532 to 15 + 6084 w 16 +■■■ 

(2) (2) (3) 

and three solutions S{ , and S{ with logarithmic terms, behaving, for small 
to, as follows : 

S[ 2) = Mw)-(S 1 -4S 2 +1QS 3 -216 59) -32w 4 -c 2 , 

Sf ' = ln(to) ■ (S 2 - 2 S 3 + 24 S 9 ) + 8 to 4 • d 2 , 

S[ 3) = 3 ln(w) 2 • (Si + 5S 2 - 2 S 3 ) 

-6 ln(to)-(2S 3 -Sf } -9S^ 2) ) -19w 4 -e 3 

where c 2 , d 2 , e 3 denote functions analytical at to = 0: c 2 = 1 + 167/96 to + 
2273/96 to 2 + d 2 = 1 + 5/2 to + 103/4 to 2 + 315/4 to 3 + ■■■, and e 3 = 1 + 
7693/456 to + 575593/11400 to 2 + 2561473/5700 to 3 + •• • 

Remark: Besides the known regular "Nickel-type" singularities (3) mentioned 
above, we remark the occurrence of the roots of the polynomial P7 of degree 28 in 
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w, which just correspond to apparent singularities [8] and the occurrence of the two 
quadratic numbers l + 3w+4w 2 = which are not of the form (3). The two quadratic 
numbers are not on the s-unit circle : \s\ = \[2 and \s\ = 1/V2. 

These new quadratic singularities that are not of the "Nickel type" (3), are 
singularities of the Fuchsian differential equation associated with x^ , but are not 
singularities of the ("physical") solution x^ itself [11] ! More generally, for arbitrary 
X^'s, the "non-Nickel" singularities [23], are singularities of the Fuchsian ODE's, but 
does not seem to be singularities [12] of the x^'s. This is a very strong motivation to 
get more informations on the next x^ 's (x^ 7 > •••)> m particular their singularities 
and the singularities of the associated [12] Fuchsian ODE's. 

The best series analysis of the x ( '"' ) ' s i which does not take into account the 
existence of the corresponding ODE, will not be able to find such singularities^ . These 
singularities are beyond any series analysis of the x^'s. 

It is probably not necessary to underline the fundamental role of the 
understanding of singularities in lattice statistical mechanics, as far as physics is 
concerned. Along this line let us mention the recent comment by B. Nickel [23] 
on our two first papers [8, 9] where he gives strong arguments that these new sets 
of singularities should actually (necessarily) correspond to pinch singularities [23] 
totally similar to the well-known Landau singularities encountered in Feynman 
diagrams [23, 24, 25, 26, 27]. Such "Landau singularities" can be seen as some frontier 
between different perturbation regimes and their implications on the physics are too 
numerous to be listed here. As far as terminology is concerned we should not call 
these pinched "non-Nickel" singularities "Landau singularities" but rather "Landau- 
like" singularities : they are singularities of a mathematical structure associated with 
the x^'s (namely the associated Fuchsian ODE) and not automatically singularities 
of the x^' s themselves. 

2.2. The Fuchsian ODE for x (4) . 

It is clear that the four particle contribution x^ is even§ in w. We thus introduce, in 
the following, the variable x = 16w 2 . The same expansion method [9] is applied with 
some symmetries and tricks specific to x' 4 ' to generate 216 terms in the variable x. 
We succeeded in obtaining the differential equation for x^ that is given in x = 16 w 2 
by 

10 d n 

E «»(*) • j-zn*) = 0) 



dx r - 

n=l 



with 



a w = -512 x e (x - 4) (1 - xf P w {x), 

a 9 = 256 {l-xfx 5 P 9 (x), a s = -384 (1 - xf x 4 P s {x), 
a 7 = 192(l-xfx 3 P 7 (x), a 6 = -96 (1 - xf x 2 P 6 {x), 
a 5 = 144 (1 -x) xP 5 (x), a 4 = -72P 4 (x), a 3 = -l08P 3 (x), 
a 2 =-54P 2 (x), ai=-27Pi(z) (10) 

ff A simple Pade analysis of x' n '' s scries expansions will not see such "non-Nickel" singularities. A 
very careful diff-Padc analysis, taking into account the Fuchsian character of these ODE's, should be 
able to see such "non-Nickel" singularities, when a straight simple diff-Padc analysis will not. 
§ This is also the case for any x^ 2 ™' (see for instance [7]). 



Beyond series expansions. 



6 



where P w (x), P 9 (x) ■ ■ ■, P\{x) are polynomials of degree respectively 17, 19, 20, 21, 
22, 23, 24, 23, 22 and 21, and are given in [10]. 



•T'-ci ri en 1 1 ft ri f v 


c-Qi n p"i 1 1 ft ri tv 


Critical OYnnnpnts in t 


p 





0,oo 


8,3,3,2,2,1,1,0,0,-1/2 


1 


1 


±1 


3,2,1,1,0,0,0,0,-1,-3/2 


3 


4 




8,7,13/2,6,5,4,3,2,1,0 





oo 


±i 


5/2,3/2, 3/2, 1/2, 1/2, 1/2, 1/2,0, -1/2, -1/2 


3 


xp, 17 roots 


sp, 68 roots 


10,8,7,6,5,4,3,2,1,0 






Table 2: Critical exponents for each regular singular point. P is the maximum 



power of the logarithmic terms for each singularity, xp is any of the 17 roots of P\q(x). 

2. 3. Seeking for the simplest Fuchsian ODE: paradox and subtleties. 

As far as finding the ODE's satisfied by the X (n) 's is concerned one can, of course, 
find several ODE's for a given ■ Actually, denoting L^™-* the differential operator 
of smallest order q (order seven for x^, ten for x*- 4 \ ••■) associated with a given 
^(n) (x^) = 0), the other ODE's can be seen to be associated with differential 
operators of order p of the form L p n ' = • Lq™' , p > q. At this step one should 
try to fight the prejudice that L p ™\ being the product of L-? and (the smallest 

order) operator Lq™\ is necessarily more complicated than these two differential 
operators: in some cases it may actually be simpler ! Among these various ODE's and 
associated differential operators, some are singled-out: the differential operator Lf 
such that one does not have any apparent singularity anymore is clearly singled-out. 
Note that the requirement to have no apparent singularities, or just a polynomial 
corresponding to apparent singularities of small degree (quadratic ...), corresponds to 
ODE's such that the number of coefficients in a series necessary to guess the ODE, 
is actually much smaller than the number required to guess the differential operator 
Lg™' corresponding to the ODE of smallest order. However, the requirement to have 
no apparent singularity is not the optimal condition to find (guess ...) the ODE with 
the minimal number of coefficients. The optimal order p for guessing the (Fuchsian) 
ODE verified by a given \^ an d its associated differential operator Lp™' 1 , is, in fact, 
an integer such that q < p < m. 

For instance, for \^ one can introduce the following differential operators 

r - t (4) - 

£14 = L 14 • £io, £13 = L}g • £io (11) 

£l2 = Li2'Ao, Al = I'll ' ^10 

where £io is the differential operator of smallest order which requires 242 coefficients 
to be guessed; £14 is a differential operator with one apparent singularity that requires 
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137 coefficients to be guessed, £13 is the differential operator which requires the 
smallest number of coefficients to be guessed namely 134 coefficients (and it has a 
quadratic apparent polynomial). £ 12 requires 143 coefficients and has a quartic 
apparent polynomial. L\\ requires 161 coefficients and it has an apparent polynomial 
of degree seven. Note that L}" is not divisible by L}° (and similarly LJ3 by L}°, etc 
...). 

Therefore, we see that there are, at least, three singled-out ODE's : the one of 
smallest order q, the one with no apparent singularity of order m, and the "optimal 
one to be guessed' which is of order p with q < p < m. 

2.4- Singled-out solutions of the Fuchsian ODE for \^ an d associated 
factorizations of differential operators. 

The order ten Fuchsian differential equation (9) has three remarkably simple algebraic 
solutions, namely the constant solution, Sq{x) — constant, and : 

_ , , 8 - 12a; + 3x 2 2-Qx + x 2 

Slix) = 8(1-303/* ' S2{X) = 2(l-s)VS - 
These solutions correspond to solutions of some differential operators of order one. 
Denote these three order one differential operators respectively in, L\ and L^. 

A truly remarkable finding)! is the following solution of the Fuchsian differential 
equation (9) : 

Ss(x) = ^ '2^1 (|| 3; x) (12) 

which is nothing but the two-particle contribution to the magnetic susceptibility, i.e., 
X^ 2 ) associated with an operator of order two, N a such that No(S 3 ) = 0. 

We do not see this result as a coincidence, but rather as the emergence of some 
"Russian-doll" structures\\ for the infinite set of the x^' s ; an d therefore as a first 
example of a structure bearing on the whole susceptibility \. 

The second solution of the order-two operator N a is simply given^f by: 

S 3 (x) = ^x 2 - 2 ^(||2;l-x) (13) 
which can also be written around x = as (^(x) analytic): 

e / n c^,^ 1 ^ k d /r(fe-i/2)r(fe + i/2)\ 

S 3 (x)= ^)-ln(,)+— V( fc -i)r(fc + l) ) 

It is obvious on the simple form (13) that this second solution is not singular at the 
ferromagnetic critical point x = 1, but is actually singular at x = 0. 
We found another solution that reads 

r, , x 4(x- 2V4 - x „„, x ,„ 

^ = x 1 + 161n ( 2 + V4^) 2 (U) 

+ 16S 1 (x)-]ng(x) -16^S 2 (x) . 2 iq(||| J 

(j Do note that we had a similar, though less spectacular, property for the differential operator of 
x' a ) : solution Si in (7) is nothing but jjjW and is thus also a solution of that operator. 
|| For instance x' 3 ' could be a solution of the differential operator associated with x^ 15 "*; but a ^ s0 of 
the differential operator associated with x' 9 '- 

1 Or in terms of the MeijerG function [28, 29] as tt/12 • MeijerG ([[], [1/2, 3/2]], [[2,0], Q],x) (see [10]). 
We thank P. Abbott for pointing out a missprint in the corresponding formula given in a previous 
version of [10]. 
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g(x) = - ((8 - 9x + 2x 2 ) + 2 (2 - x) ■ y/(l - x)(4 - xfj (15) 

The differential equation corresponding to Si(x) has also So, Si(x) and S 2 (x) as 
solutions. 

Denoting by Ciq the differential operator corresponding to (9), one may search 
for simple solution of its adjoint do- We found a solution corresponding to an 
operator of order two (denoted N%) which is a combination of elliptic integrals with 
polynomials [10] of degree 21 

where : 

/f(a;)= 2^(^:1; a:), %) = 2 Fi(-^;l ;i ). 

All these solutions of (9) can be used to build the complete factorization of the 
differential operator do that we write as 

C w = N$ ■ M 2 ■ G(L) = Mi • N 9 ■ G(L) = M t ■ L 24 ■ G(N) 

where G(L) is a shorthand notation of a differential operator of order four, factorizable 
in four order one operators)} and where G(N) is an operator of order five, factorizable 
in one operator of order two and three operators of order one. G(N) has 24 different 
factorizations involving eight differential operators of order two, and 24 order one 
operators with twelve appearing in G(L). All these 36 factorizations (16) are given in 
[10]. 

From the 24 factorizations of the order five differential operator G(N), six are 
divisible by the differential operator N . One such factorization reads G(N) = 
Li 3 ■ Lyj ■ Ln ■ N a implying the direct sum decomposition : 

do = C 8 ®N a , C 8 = M 2 -G(L) (16) 



2.5. More factorizations for the differential operator of ■ 

Let us come back to x^. One can changeff S3 into S^ new ^ — S3 — 16 5*9. This 
amounts to changing the w 9 coefficient in S3 from 7479 to 7463 (and the following 
coefficients consequently). Actually B. Nickel foundj that the new solution S^ new ^ is 
actually a solution of a fourth order differential equation. 

Focusing on this differential operator one can^f find the general solution of this 
order four differential operator, which we call N4, which reads: 

Sol = ot\ ■ Si + a 2 ■ S 2 + S2 • f 2 ■ dw (17) 

U The order four operator G(L) has six different factorizations involving thirteen differential operators 
of order one. 

ff In some of our calculations the linear combination 53 — 12 Sg which amounts to changing the u> 9 
coefficient in S3 from 7479 to 7467 has also been used. 
X Private communication. 

K Using for instance the dsolve and DFactor command of DEtools [30]. See also [31, 32, 33]. 



Beyond series expansions. 



9 



where Sol 2 denotes the general solution of an order two differential operator [11] Z 2 : 
_ d Ad B 
2 dw 2 P dw P ' 
P = w ■ (1 - 4 wf (1 + 4 w) (1 + 3 w + Aw 2 ) (1 - w) (1 + 2 w) ■ Pi 

Pi = 96 w 4 + 104 w 3 - 18 w 2 - 3 w + 1 

A = (1-Aw) (49152 w 10 + 72704 w 9 - 48384 w 8 - 94272 w 7 

- 40368 w 6 - 4488 w 5 + 1080 w 4 - 108 w 3 - 111 w 2 - 6 w + l) 

B = 98304 w w + 98304 w 9 - 239616 w s - 307200 w 7 - 71552 w 6 
- 4224 w 5 - 3192 w 4 - 1520 w 3 - 276 w 2 + 48 io + 4. 

Let us make a first remark: the degree four polynomial P\ corresponds to 
apparent singularities of the order two differential operator Z 2 - This differential 
operator appears when solving the order four differential operator N4 which has seven 
apparent singularities roots of a polynomial of degree seven. One thus sees that the 
degree 28 polynomial (denoted P7 in [8], see equation (9) in [8]), corresponding to 
28 apparent singularities, has been replaced, as far as ./V4 is concerned, by a degree 
seven polynomial or, as far as Z 2 is concerned, by the degree four polynomial Pi. The 
apparent singularities of these various differential operators are clearly highly volatile. 

The linear combination (17) of four simple solutions is in complete agreement 
with the (Dfactor/DEtools) factorization of this order four differential operator 7V4: 

N A = N 3 ®L U Li (Si) = (18) 

iV 3 - Z 2 -N u JVi(S 2 ) = 0. 

From the existence of this order four differential operator one immediately deduces 
an order three differential operator L3, such that the order seven differential operator 
I/7 factorizes as follows: 

L 7 = L 6 £1, with: L 6 = L 3 ■ N 3 = L 3 ■ Z 2 ■ N x (19) 

Let us note that the four solutions (17) are totally decoupled from the linear 
combination Sg + ji ■ S\ such that Lq(Sq + n ■ S\) = 0. 

Let us also mention that at the singular points w = 1, w = — 1/2, and at the two 
quadratic roots w\, w 2 of l + 3w + 4w 2 = 0, the solution carrying a logarithmic term is 
in fact a solution of Z 2 -N\. Therefore, the three solutions of the differential operator 
L 3 ■ Z 2 ■ N\, emerging from L 3 , are analytical at the non physical singular points 
w = 1, w = —1/2, and at the quadratic roots of 1 + 3w + 4w 2 = 0. The solutions of 
the differential operator L 3 , itself, are of "elliptic integral type" [11]. At the singular 
point w = 1/4, one also notes that the differential operator Z 2 ■ N\ is responsible of 
the (1 — 4u>) _1 behavior. One will then expect the "ferromagnetic constant" 1$ to 
be localized in the blocks of the connection matrix corresponding to the solutions of 
the order three differential operator Z 2 ■ Ni at the point w = 1/4. 

Note that finding such complete factorization for the differential operator L-j (or 
£10 for x^) is °f great help to get the differential Galois group (which identifies here 
with the monodromy group), and, more precisely, when computing the connection 
matrices necessary to write the monodromy matrices written in the same (Frobenius 
series solutions) basis. 
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3. Differential Galois group, singular behaviour and asymptotics 

3.1. Differential Galois group. 

From this factorization (18) of Lj one can deduce, immediately, the form^f in some 
"well-suited global basis" of solutions, of the 7x7 matrices representing the differential 
Galois§ group [11, 35, 36, 37] of L-j, namely: 



al 























bl 




















hi 


gi 


g2 














h2 


gs 


gi 














HI 


H2 


H3 


Gl 


G2 


G3 





H4 


H5 


H6 


G4 


G5 


G6 





HI 


H8 


H9 


Gl 


G8 


G9 



(20) 



where the 2x2 and 3x3 block matrices Qal(Z 2 ) and Qal(L 3 ) 



Qal{Z 2 ) 



gi g2 
g3 g4 





' Gl 


G2 


G3 


Qal{L 3 ) = 


G4 


G5 


G6 




Gl 


G8 


G9 



correspond respectively to the differential Galois groups of the differential operators 
Z 2 and L 3 . The 3x3 matrix Qal(N 3 ), corresponding to the differential Galois group 
of N 3 = Z 2 N\, and the 3x3 block- matrix corresponding to the fact that we have 
a semi-direct product of the differential Galois group of L 3 and N 3 in Lq = L 3 ■ N 3 , 
read: 

"WOOl f HI H2 H3 



Gal(N 3 ) = 



hi gl g2 
h2 g3 g4 



H(N 3 ;L 3 ) = 



H4 H5 H6 
HI H8 H9 



The semi-direct product structure is also clear on the 3x3 matrix Qal(N 3 ). The 
3x3 matrix block H(N 3 ;L 3 ) (rcsp. the entries hi and h 2 in (20)) can, in principle 
be calculated exactly from|| the differential Galois group of N 3 , and the differential 
Galois group of L 3 (resp. the differential Galois group of L 2 and the differential 
Galois group of N\). The calculation of these off-diagonal blocks, associated with the 
fact that we have a semi-direct product instead of a product (a product of differential 
operators instead of a direct sum of differential operators), relies on the fact that one 
knowns how to actually calculate the morphisms between differential operators (see 
for instance P. Berman and M. Singer [39] or P. Berman [40]). ft is worth noting 
that such calculations are totally effective. We will, not however, give more details of 

% One has to be careful with the definition one takes for these matrices, which will be associated 
with the choice of a vector set of solutions, as a line vector or as a column vector. Different choices 
change these matrices into their transpose. 

§ Since our ODE is a Fuchsian ODE the differential Galois group identifies with the monodromy group 
generated by the monodromy matrices written in the same basis: there are no Stokes matrices [34] 
associated with irregular singular points. 

| In mathematical words the solutions of L3 are in the Picard-Vessiot extension of Lg. 
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these later calculations [39, 40, 38]: they are only necessary if one wishes to perform 
exhaustively the complete analysis of the differential Galois group of L7: we just 
perform, here, a "warm-up", only sketching the analysis of the differential Galois 
group of L7, which reduces (up to a product by the field of complex numbers C, 
because L7 is the direct sum of L§ and L\) to the differential Galois group of Lq, 
which is, thus, just seen as a semi-direct product of the differential Galois group of L3, 
of the differential Galois group of Li and of the differential Galois group of N\ (namely 
the field of complex numbers C). At first sight the differential Galois group of L3 can 
be either the group SX(3, C), or the group! PSL(2, C), or the semi-direct product 
of C and of the permutation group S 3 of three elements. Some calculations that 
will not be detailed heref, show that the differential Galois group of L3 is actually 
SL(3, C) because the Wronskian of L 3 is rational and that the differential Galois 
group of L 2 is SL(2, C) also because its Wronskian is rational. The differential 
Galois group of L7 thus reduces (up to a product by C) to semi-direct products 
of SL(3, C), SL(2, C) and C. This is a quite general description of the differential 
Galois group, that may satisfy mathematicians, but, as far as physics is concerned, 
we need to actually calculate specific elements of the differential Galois group, in 
particular the "non-local" connection matrices, in order to have some understanding 
of the analytical properties of the solutions [11] (dominant and sub-dominant singular 
behaviours, asymptotics, see section (3.3) below). As far as explicit calculations 
are concerned, we have, however, learned, from this quite general analysis, that a 
well-suited "global" basis for writing explicitely any kind of " non-local" connection 
matrices, should be : Si, S2, s[ 2 ' new \ S^ new ^ (corresponding to the solutions of (17), 
on the N4, side) together with the further solutions S9, S[ 3 ^ coming from the 

right division of L7 by N4. 

In particular it is shown in [11], and this will be sketched in the next section 
that the connection matrix between the singularity points and 1/4 (matching the 
well-suited local series-basis near w = and the well-suited local series-basis near 
w = 1/4) is a matrix of the form (20) where the entries are expressions in terms of 
ln(2), ln(6), V3, tt 2 , 1/tt 2 , ... and a constant 1^ introduced in equation (7.12) of [1]: 

I I dyidwdK, ((y?-1)(i-l)) 1/2 - y2 = 
= .000814462565662504439391217128562721997861158118508 • • • 



(2/1 + 2/2) (2/2 + 2/3) (yi +V2 + yz) ' 

This constant can actually be written in term of the Clausen function C7 2 : 

J 3 + - i"(y +2-3V3.C*4)) (2D 
where CI2 denotes the Clausen function : 

n=l 

Similarly one can consider the (Frobenius series) solutions of the differential 
operator associated with around x = and around the ferromagnetic and 

J Such a PSL(2, C) situation would correspond to building differential equations such that one of 
their solutions is actually the product of two solutions of a second order differentiable equation. One 
can easily build such examples using the diffeq * diffeq command in Gfun [41]. 
f A publication with Jacques- Arthur Weil on these questions will follow [42], 
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antifcrromagnetic critical point x = 1 respectively. Again the corresponding 
connection matrix (matching the solutions around the singularity points x — and 
the ones around the singularity point x = 1) have entries which are expressions in 
terms of 7r 2 , rational numbers but also of constants like constant 1^ introduced in [1] 
which can actually be written in term of the Ricmann zeta function, as follows : 

The derivation of these two results (21), (22) has never been published^ but these 
results appeared in a conference proceedings [43]. We have actually checked that 1$ 
and we got in our calculations of connection matrices displayed in [12] (and in the 
next section) as floating numbers with respectively 421 digits and 431 digits accuracy, 
are actually in agreement with the previous two formula. These two results (21), (22) 
provide a clear answer to a question that will be addressed in the next section of how 
"complicated and transcendental" some of our constants occurring in the entries of 
the connection matrices can be. These two remarkable exact formulas (21), (22) 
are not totally surprising when one recalls the deep link between zeta functions, 
polylogarithms and hypergeometric series [44, 45, 46, 47]. Along this line, and keeping 
in mind that we see all our Ising susceptibility calculations as a "laboratory" for other 
more general problems (Feynman diagrams, ...), we should also recall the various 
papers of D. J. Broadhurst [48, 49, 50, 51] where C^(§) and £(3) actually occur in 
a Feynman-diagram-hypergeometric-polylogarithm-zeta framework (see for instance 
equation (163) in [48]). 



3.2. Connection matrices and monodromy matrices. 



Connecting various sets of Frobenius series-solutions well-suited to the various sets 
of regular singular points amounts to solving a linear system of 36 unknowns (the 
entries of the connection matrix). We have obtained these entries in floating point 
form with a very large number of digits (400, 800, 1000, ...). We have, then, been able 
to actually "recognize" these entries obtained in floating form with a large number 
of digits [11]. As an axample the connection matrix for the order six differential 
operator Lq matching the Frobenius series-solutions around w — and the ones 
around w = 1/4, namely C(0, 1/4) reads: 



C(0, 1/4) 



(23) 



1 














■ 


1 





9V3 

64tt 














3tt\/3 
32 














5 


3 Z J 3 


3^ 
64tt 








1 

16tt 2 


5 
4 


371-^3 
32 


45^3 
256tt 





1 

32 





27T 2 

3 


157rV3 
64 


225\/3 37rV3 
1024tt 64 


64 









29 
16 

Not surprisinglyfl a lot of 7r's pop out in the entries of these connection matrices. We 
will keep track of the 7r's occurring in the entries of connection matrices through the 

% We thank C. A. Tracy for pointing out the existence of these two results (21), (22) and reference [43] . 
ft One can expect the entries of the connection matrices to be evaluations of (generalizations of) 
hypergeometric functions, or solutions of Fuchsian differential equations. 
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introduction of the variable a — 2 in. 

From the local monodromy matrix Loc(ft) , expressed in the w = 1/4 well-suited 
local series-basis, and from the connection matrix (23), the monodromy matrix around 
w = 1/4, expressed in terms of the (w = 0)-well suited basis reads: 

A 

24a 4 -M w=0 (l/4)(a, ft) = 



B C j 

C(0, 1/4) ■ Loc(fi) ■ (7(0, 1/4)- 1 



(24) 



where 



A 
B 



and [ C ] read respectively: 



-24 a 4 
-48 a 4 


-48 pi 
12a 2 p 3 
(87 + 8 a 2 )a 4 





24 a 4 


32ftp 2 
4(75-4a 2 )a 2 ft 






-144 a 2 ft 
24 a 4 
48 ft (9 a 2 +80 ft) 

-300a 2 ft 
3(4a 2 -75)a 2 ft 



with pi = 5 a 4 + 8 ft 2 + 8 ft 2 a 2 , p 2 = 4 ft a 2 - 75 ft - 15 a 2 and p 3 = 5 a 2 + 4 ft + 
4 ft a 2 , and: 





' 24 a 4 


-384a 2 ft 


1536 ft 2 


c = 





24 a 4 


-192a 2 ft 










24 a 4 



Note that the transcendental constant I 3 has disappeared in the final exact expression 
of (24) which actually depends only on a and ft. This (a, ft) way of writing the 
monodromy matrix (24) enables to get straightforwardly the N— th power of (24): 

M lo=0 (l/4)(a ) Sl) N = M tu=0 (l/4)(a, TV • ft) 

Our "global" (800 digits, 1500 terms in the series, ...) calculations yield 
quite involved exact connection matrices. With such large, and involved, computer 
calculations there is always a risk of a subtle mistake or misprint. At this step, and in 
order to be "even more confident" in our results, note that the monodromy matrices 
must satisfy one matrix relation which will be an extremely severe non-trivial check 
on the validity of these eight matrices M,, or more precisely their (a, ft) extensions. 
Actually it is known (see for instance Proposition 2.1.5 in [37]), that the monodromy 
group^ of a linear differential equation (with r regular singular points) is generated 
by a set of matrices 71, 72, • • • , 7 r that satisfy 71 • 72 • • • 7r = Id, where Id denotes 
the identity matrix. The constraint that "some" product of all these matrices should 
be equal to the identity matrix, looks quite simple, but is, in fact, "undermined" by 
subtleties of complex analysis on how connection matrices between non neighboring 
singular points should be computed. The most general element of the monodromy 
group reads: 

(25) 



1V1 P(l) 1V1 P(2) IVI P(3) 



IVI P(4) IVI P(5) IVI P(6) IV1 P(7) IVI P(8) 



K Which identifies in our Fuchsian case to the differential Galois group; we have regular singular 
points, and no irregular points with their associated Stokes matrices [34]. 
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where P denotes are arbitrary permutation of the eight first integers. In other words, 
one of the products (25) must be equal to the identity matrix for some set of rij's 
and for some permutation P. Let us introduce the following choice of ordering of the 
eight singularities, namely oo, 1, 1/4, w\, —1/2, —1/4, 0, w 2 (w\ = (3 + i-\/7)/8 
and W2 = w* are the two quadratic number roots of 1 + 3w + Aw 2 = 0), the first 
monodromy matrix M\ is, thus, the monodromy matrix at infinity M.(co), the second 
monodromy M2 matrix being the monodromy matrix at w — 1, M(l), ... This is 
actually the particular choice of ordering of the eight singularities, such that a product 
like (25) is equal to the identity matrix^f: 

Mi • M 2 • M 3 • M 4 • Ma ■ M e ■ M 7 ■ M 8 = Id (26) 
= M(oo)-M(l)-M{l/4:)- M(w!) 

x M(-l/2) • M(-l/4) • M(0) -M(w 2 ) 

It is important to note that these relations (26) are not verified by the (a, 17) extension 
(like (24)) of the monodromy matrices Mj. If one consider relation (26) for the (a, 17) 
extensions of the Mi's, one will find that (26) is satisfied only when a is equal to 
17, but (of coursefl) this a = 17 matrix identity is verified for any value of 17, not 
necessarily equal to 2 in. 

A last remark is the following: right now, we have considered all the matrices 
(connection and therefore monodromy matrices expressed in a unique non-local basis) 
with respect to the (w = 0) well-suited basis of solutions. This is motivated by the 
physical solution x*- 3 ^) which is known as a series around w = 0. In fact, we can switch 
to another w = w well-suited basis of solutions. This amounts to considering the 
connection C(w, w s ) = C (0, w) ■ C(0, w s ). For instance, we have actually performed 
the same calculations for the (w = 1/4) basis of series solutions. We have calculated 
all the connection matrices from the (w = 1/4) basis to the other singular point 
basis series solutions, and deduced the exact expressions of the monodromy matrices 
now expressed in the same (w — 1/4) basis of series solutions. It is worth noting 
that we get, this time, for the monodromy M u)=1 / 4 (w s ) around the singular point 
w s and expressed in the (w = 1/4) basis, a matrix whose entries depend rationally 
on a — 2 in, 17 = 2 in, but, this time, also ( except for the monodromy matrix for 
w = \) on the "ferromagnetic constant" 1% . One verifies that the product of these 
(a, f2, /.^-dependant monodromy matrices in the same order as (26), is actually equal 
to the identity matrix when a = 17, the matrix identity being valid for any value of 
a = 17 (equal or not to 2 in), and for any value of l£ (equal or not to its actual 
value) . 

We have similar results for the monodromy matrices around singular point w s , 
expressed in the (w = 00) basis, but, now, the monodromy matrices M w=ao (w s ) 
depend on a = 2 in, 17 = 2 i n, and, this time, on some (as yet unrecognized) 
constants j/41 and X42 (which can be found in the Appendix C in [11]). Again, the 
product of these monodromy matrices in the same order as (26), is actually equal to 
the identity matrix when a = CI, the matricial identity being valid for any value of 
a = 17 (equal or not to 2 in) and for any values of j/41 and X42 . 

If Of course, from this relation, one also has seven other relations deduced by cyclic permutations. 
(J A matrix identity like (26) yields a set of polynomial (with integer coefficients) relations on 
Q = 2 in. The number n being transcendental it is not the solution of a polynomial with integer 
coefficients. These polynomial relations have, thus, to be polynomial identities valid for any Q. 
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3.3. From the differential Galois group to the singularity behaviour of solutions and 
other asymptotics \ studies. 

The knowledge of these connection matrices gives, in particular, the decomposition of 
the ' s on some series well-suited to the analysis of the neighborhood of the critical 
points and which, consequently, have a well-defined critical behaviour near the critical 
points. One can, then, immediately deduce the singularity behaviour of the x^' s - 

Let us focus, for instance, on the critical behaviour of x' 3 ^ near the ferromagnetic 
critical point w — 1/4. Denoting z = 1 — Aw, the singular part of the "physical" 
solution Sg, associated with x^ 3 \ reads : 

j+ 

X^ (singular) — 8 • Sg(singular) ~ + ln(z) • Singi 



64 7T 



(ln(z) 2 - 6 ln(2) m(z)) • Sing 2 



where 



Sin gi = 23/3 + 227^/240 -2047z 2 /1680 - 88949 z 3 /92160 +••• 

Sing 2 = 1 - z/8 + 3 z 2 /16 + 29 z 3 /512 + ■ ■ ■ 

and where l£ is actually the transcendental constant (21), namely 1$ = 
0.000814462 • • • These results agree with previous results of B. Nickel (see for instance 
the footnote of page 3904 in [6]), but the correction terms are new§, in particular the 
terms in ln(2)/(647r 2 ). In the r = (1/s — s)/2 variable introduced in [15, 16] this 
reads 

where the dots denote less divergent terms. Near the antiferromagnetic critical point 
w = — 1 /4 we get : 

X {3] 1 (singular) ~ — \ ■ ln 2 (l + 4 w) + ■ ■ ■ (27) 

3271" 

From these results, and using some identities likefj 
ln 2 (l -x)= J2 ~ ■ (^( n ) + l)-x n 

where 7 denotes the Eulcr Gamma function, and ip(n) is the logarithmic derivative of 
the r function, and recalling, that, when n — > 00, one has the asymptotic behaviour 
ip(n) — > ln(n), one can find the following asymptotic behavior for the expansion 
coefficients c(n) of x^/8/w 9 near w — : 



c(n) ~ -i- 



/+ 1 ln(n) t b x 



4 8 • 4" w 2 32tt 2 n n 

'bi 1 ln(n) 



+ (-!)"•(-- 
V n 



+ 



n 16ir 2 n 

with 61 ~ 0.0037256 and b 2 — 0.0211719 . This parity behavior explains why we had 
so much difficulty to actually find the "true" asymptotic behavior of these coefficients 

§ All these results have also been found by B. Nickel (private communication). 

J The asymptotic expansion of the coefficients of Taylor series can be obtained using various packages 
availables at http://algol.inria.fr/librarics/softwarc.html. For instance one can use the command 
"equivalent" in gfun [41] (see [52] for more details. For more details on calculations of asymptotics 
see [53, 54, 55]. 
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after obtaining around 500 coefficients (see equations (3.6) in [9]: one has 13.5 x 4 
for N even and 11 x A N for N odd, instead of the "true" (4 7 • 7+) x 4 N ~ 13.4 x 4 
asymptotic behavior). 



5.^. Generalization to x^ 

All the previous calculations can be performed, mutatis mutandis, for x^- Fot 
instance from the singular behaviour near x = 1 : 

X^(sm^ar) ~ + — j-^ • ln 3 (l - x) + ■ ■ • (28) 

where J 4 ~ is given by (22) ( J 4 _ = 0.00002544851106586 • • • ), one obtains the following 
asymptotic form for the series coefficients)! c(n) of x^ near x = : 

1 ln 2 (n) ln(n) a 3 

c n ~ J 4 " - — — • — ^ + a 2 ■ + -i + • • • 29 

128 7r d n n n 

where a 2 ~ 0.0012515 and a 3 ~ -0.0021928. 



3.5. Non generic character of 

As we remarked in [9] the critical behavior of x^ near the ferromagnetic point 
w = 1/4, does not corresponds to the "dominant" exponent one obtains in the set of 
indicial exponents of the indicial equation at w = 1/4. The same situation occurs 
exactly for x^ as wm be seen below. The singular behavior of the solutions of the 
differential equation can be read easily from Table 2. Near x = 4, they are i 13 / 2 , 
where t = A-x. Near x = oo, they are t~ 1/2 log fe (t), (with k = 0, 1), t 1/2 ln fc (t) (with 
k = 0, 1,2,3). t 3/2 ln fc (t) (with fc = 0, 1) and £ 5 / 2 , where t = 1/x. Near x = 1, they 
are t~ 3/2 , i -1 , ln fe (t) (with fc = 1,2,3) and f ln(i), where t = 1 - x. Note that these 
behaviors take place for the general solution of the differential equation (9). 

As far as the physical solution is concerned, the dominant singular behavior at x = 
4 (namely i 13 / 2 ), and x = oo (namely i -1 / 2 ln(i)), are present in the physical solution 
X^ confirming [7, 15]. At the ferromagnetic point (x = 1 which also corresponds to 
the anti-ferromagnetic point for x^ ): with the dominant behavior i -3 / 2 , the growth of 
the coefficients would be (3/2)^ /AM <~ >/N. Since this is not the case, the coefficients 
of the series for large values of N behaving like CV ~ 0.2544 • 10~ 4 , this t~ 3 / 2 
singular behavior does not contribute to the physical solution. Only the subdominant 
singular behavior is actually present. 

As far as singularity analysis in lattice statistical mechanics and, more specifically, 
as far as the analysis of singular behavior near a physical singularity, like the 
ferromagnetic singularity x — 1, is concerned, we probably have a prejudice which 
consists in believing that "physical solutions" will correspond to generic solutions and, 
thus, will take into account the strongest possible singularity if a set of exponents is 
available (from the indicial equations). In fact the "physical solution" x'" 4 '' chooses to 
"avoid" this dominant singular behavior, and will, near the ferromagnetic singularity 
x = 1, correspond to "sub-dominant" singular behavior. 

As far as asymptotic calculations are concerned, there are many tools and 
methods [53, 54, 55] which give (after some work) asymptotic results for the "generic" 
solutions of an ODE (or of a recursion). In the "sneaky" situation where physics 

D More generally terms like ln l (l — x) give coefficients growing asymptotically like 
(-ir(ln(n))(-D/n. 
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does not correspond to the "generic" solutions associated with the strongest possible 
singularities, all the asymptotic calculations become much more difficult. In fact such 
calculations really correspond to this difficult class of "global" problems of finding a 
"non-local connection matrix" . 



3.6. + versus x at scaling. 

Thus far we have sketched some mathematical aspects of the solutions to the Fuchsian 
differential equations for x^ and X^- However, the physics implications of the 
solutions we have obtained call for some remarks near the physical critical points. 
Taking, as an example, the ferromagnetic singularity for x' 3 \ the sum of the first two 
n-particle terms behave near r = (1/s — s)/2 ~ as: 



,<., + ^i±i_^ + ^.g (30) 



11 3 , 1 /, 2/ „, 23 , 14 \ 

H h -it — n 2 2 n(2 H M 

48 8 3 4tt 2 V 12 144/ 

The exact susceptibility, as reported in [16], yields for the normalized susceptibility x 
near r = (1/s - s)/2 ~ 0: 

X (1- S 4 )V4' X (l+T 2 ) 1 ^ X [6L > 



— 1/4 oc oo 

; Cl r- 2 F + (r) + ^ ^ ' T? ln ^ T ) 

p— g— p 2 



where c\ = 1.000815260- • • is given with some 50 digits in [15]. F + (t) and b+' 9 ^ arc 
given in [15]. The constants 1 + 1$ and c\ verify 1 + 1$ + l£ = c\ with nine digits, 
I§, corresponding to x' 5 \ IS a constant given in [1] (and with some 30 digits in [6]). 
Thus, and as suggested in [1], the partial sums of the x^ should converge rapidly 
to the full X- Furthermore, adding x^ term has resulted in a series expansion that 
reproduces the first 24 terms of x to be compared with only eight first terms for x^ 
series. 

However, equation (31) shows a r -1 / 4 divergence as an overall factor to the 
logarithmic singularities which is absent in (30). The same situation occurs for the 
low temperature regime when we compare the first two n-particle terms (x*' 2 *' and x*" 4 "*) 
with the full x at scaling ff . This observation raises several profound issues: how does 
the logarithmic terms in the entire sum add up to make the r -1 / 4 divergence be 
factored out ? If one assumes that the other x( 2n+1 ) terms share the same singularity 
structure as x 3 , hi particular the occurrence (in variable r or s) of only integer critical 
exponents at the ferromagnetic critical point, the t -1 / 4 divergence, as an overall factor, 
implies the following correspondence : 

oo N(n) oo oo 

£ £ a n , m ■ S n , m (r) ln m (r) r^ 4 . £ £ b ^ T » ln'(r) 

n— 1 m=0 p—0q—p 2 

where 5„^„(r) is analytical at r = 0, the a„ jm 's are numerical coefficients and 
N(n) is the maximum power of logarithmic terms occurring in the solution around 

ft For the leading amplitude, \^ an d g ive 1/127T + I± ~ 1.0009593 ■ • ■ /12n which is very close 
to 1.0009603 ■ ■ ■ /12n for the full £ [6]. 
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the ferromagnetic point of the differential equation of \ 2n+1 . This correspondence 
requires probably a very particular structure in the successive differential equations. 
Note that the discrepancy phenomenon we discuss here may be more widespread than 
that observed here. The challenging problem one faces here is to link linear and non 
linear descriptions of a physical problem, namely the description in terms of an infinite 
number of holonomic (linear) expressions for a physical quantity of a non linear nature. 
The difficulty to link holonomic versus non-linear descriptions of physical problems 
is typically the kind of problems one faces with the Feynman diagram approach of 
particle physics. 

4. Conclusion 

Most of the results we have displayed in this paper are, in large part, mathematical 
"subtleties" associated with the structure of the differential operators associated with 
the X^' s an d also surprising results like the fact that the Fuchsian ODE's of 
(resp. x^) have simple rational and algebraic solutions, together with previously 
encountered solutions x^ (resp. x^ 2 -*). 

In a miscellaneous list of these mathematical "subtleties" let us recall the 
occurrence of a large set of apparent singularities for the ODE of smallest order (or 
the fact that the simplest, and easiest to be guessed, ODE is not the one of smallest 
order), the existence of "Landau- like" pinched singularities, (that are not of the "Nickel 
form" (3)), that are singularities of the ODE but not singularities of the "physical" 
solutions x^ n \ the fact that the "physical" solutions are not "generic" solutions 
but "mathematically singled-out" ones (subdominant critical behavior§, no singular 
behavior for the previous "Landau- like" pinched singularities), ... 

The fact that x*- 2 -* is a solution of the differential operator associated with \^ is 
probably the most important result, since it is a strong indication of some "Russian- 
doll-like" mathematical structure valid for the infinite set of differential operators 
associated with all the infinite sequence of X^' s - 

It is worth noting that almost all these mathematical structures, or singled-out 
properties, are far from being specific of the two-dimensional Ising model : they also 
occur on many problems of lattice statistical mechanics or, evenff, as A. J. Guttmann 
et al saw it recently, on enumerative combinatorics problems like, for instance, the 
generating function of the three-choice polygon [56] . 

With this specific study of the susceptibility of the Ising model, we think that, 
in the future, we should contemplate a "collision" between various mathematical 
domains and structures, namely the theory of isomonodromic deformations (Painleve 
equations and their generalizations with unit circle natural boundaries like Chazy 
III equations, see [57, 58], or Gamier systems! [59]), the theory of holonomic 
systems (Fuchsian equations and in particular the so-called rigid local systems [60] 
and their geometrical interpretations [61]), the theory of modular forms [62, 63] 
(Eiscnstein series [64, 65], Heegner numbers [66], complex multiplication for elliptic 

§ At the ferromagnetic critical point w = 1/4, the behavior (1 — 4mj)~ 3 < /2 corresponding to the largest 
critical exponent for the ODE is actually absent in the physical solution x ■ 

ff The wronskian of the corresponding differential equation is also rational, the associated differential 
operator factorizes in a way totally similar to the one described in sections (2.4), (2.5), large 
polynomial corresponding to apparent singularities also occur, ... 

X The occurence of hypergeometric functions of several complex variables in the Gamier systems is 
worth noting (see page 10 and page 111 paragraph 31 in [59]). 



Beyond series expansions. 



19 



curves [66], modular equations [68, 69, 70] ...), and various ideas related to the 
analysis of Feynman diagrams (Landau singularities, generalizations of hypergcometric 
functions to several complex variables, Appel functions, Kampe de Feriet, Lauricella- 
like functions, polylogarithms [48, 49, 50, 51], Riemann zeta functions, multiple zeta 
values, ...). We also think that such "collisions" of concepts and structures are not 
a consequence of the free-fermion character of the Ising model, and that similar 
convergence should also be encountered on more complicated Yang-Baxter integrable 
models§, the two-dimensional Ising model first popping out as a consequence of its 
simplicity. An open question is what could remain of these related structures for 
less Yang-Baxter integrable problems such as the calculations of Feynman diagrams 
? Apparently a lot seems to remain [48, 49, 50, 51]. Another open problem which 
requires to be addressed is the (apparent ?) scaling discrepancy, we sketched in section 
(3.6), between the holonomic approach and more Painleve oriented approach. 
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